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SUMMARY 

The Ceramics Analysis and Reliability Evaluation of Structures (CARES) integrated ^ign program 
on statistical fast fracture reliability of monolithic ceramic components is enhanced to include the us 
neutral data base, two-dimensional modeling and variable problem size. The data base allows for the effi- 
cient transfer of element stresses, temperatures and volumes/areas from the finite element output to > the 
reliability analysis program. Elements are divided to insure a direct correspondence between the su 
elements and the Gaussian integration points. Twodimensional modeling is accomplished by assessing 
the volume flaw reliability with shell elements. To demonstrate the improvements in the algoritnm, 
example^ro^en^are selected Iron, . round-, obi* conducted by WELFEP (WEokeet L,nk Induce prob- 
ability prediction by Finite Element Postprocessors). 


INTRODUCTION 

Advanced gas turbine technologies require materials that exhibit low density, good oxidation and 
corrosion resistance, high specific strength and stiffness all at high temperature es Se l^ted u^r 

ceramics exhibit these desirable material properties. Thus, improvements in fuel efficiency are anticipa 
ed due to increased engine temperatures and pressures, which in turn generate ^eP^eran . 

Progress in utilizing advanced gas turbines has also supported cogeneration technologies. The use 
monofithic ceramics as structural materials has been limited by their inherent brittleness, susceptibility to 
damage from thermal shock, and a large variation in strength which is reflected in diminished component 

reliabilities. 

Thus, analyses of components fabricated from ceramic materials require a departure from the usual 
deterministic design philosophy (i.e., the factor of safety approach) prevalent m analyzing metallic struc- 
turluomponents For this reason, the integrated design program CARES (Ceramics Analysis and Relia- 
bility Evaluation of Structures) has been developed at NASA Lewis Research Center (Nemeth et al., 
1990) In the CARES program, a probabilistic approach is adopted that allows the design engineer o 
account for variability in strength and diminished load capacity with increasing component size (the so- 
caUedsize effect) . The goal is to predict the time-independent reliability of a structural component sub- 
jected to multiaxial load conditions. Volume and/or surface flaws can be taken into account. Wei u 


*NASA Resident Research Associate at Lewis Research Center. 

“National Research Council— NASA Research Associate at Lewis Research Center. 



material strength parameters, the Batdorf crack density coefficient, and other related statistical quan- 
tities can be estimated from bend bar or uniaxial tensile specimen fracture strength data. 

The structural response of isotropic components subjected to thermomechanical loads is usually deter- 
mined by commercial finite element programs. These finite element programs must be coupled to 
CARES for evaluation of the structural reliability for a component. The first version of CARES has 
proven its applicability in the design of gas turbine rotors, uncooled combustors, engine exhaust valves, 
laser windows on test rigs, and spacecraft activation valves. 

Currently, a next generation algorithm, CARES 3.0, has been developed that allows the design engi- 
neer access to additional finite element programs through the use of a formatted neutral data base. The 
algorithm is split into two separate modules: (1) the finite element data interface program, and (2) the 
reliability evaluation algorithm. In the first module, results from finite element analysis (available in the 
form of standard output files and/or plot data files) are interpreted. The finite element data is written to 
a neutral data base containing element or subelement stresses, volumes (or surface areas for shell 
elements), and temperatures. At this time, CARES 3.0 is coupled to the MSC/NASTRAN and ABAQUS 
general purpose finite element packages, and interfaces with ANSYS and MARC are being planned. 

In the second module, CARES 3.0 reads this data base for the reliability analysis. In the reliability 
evaluation algorithm a modular structure is used. The main program subsequently calls special sub- 
routines for reading the control file, calculating and/or storing the material parameters, printing the 
input data, reading the neutral data base, calculating risk of rupture intensities (a local measure of relia- 
bility) and probability of failure, creating a PATRAN element file, and printing the results. This modu- 
lar structure will make the implementation of future changes and additions, such as more sophisticated 
failure models, easier without undertaking extensive programming efforts. 

Furthermore, a subelement technique has been implemented. In a previous version, CARES2, the ele- 
ment stress tensor was obtained by averaging the stress state at each node over the entire element. The 
nodal stress tensor was computed by extrapolating the stresses from the integration points, and as a con- 
sequence, the stress state could only be approximated. This results in a lack of accuracy in problems 
where high stress gradients exist over the element. 

In the subelement technique the reliability function is applied at each Gaussian integration point 
instead of using averaging techniques. This approach is now feasible since all commercial analysis pro- 
grams make integration point stresses and temperatures available, although the current release of MSC / 
NASTRAN (version 66) has an error retrieving these stresses. Consequently we had to interpolate nodal 
stresses when MSC/NASTRAN was used. Thus, considerable improvements in the accuracy of reliability 
analyses have been realized. To demonstrate the improvements in the algorithm, example problems are 
selected from a round-robin conducted by WELFEP (WEakest Link failure probability prediction by 
Finite Element Postprocessors). This group includes researchers from Europe and the United States who 
are involved m the implementation of postprocessors for reliability prediction. The goals of this study are 
to compare and verify postprocessors and to formulate guidelines for their implementation and use. Ini- 
tial comparisons will be based on the results of the round-robin calculation which consists of three test 
cases. 


RELIABILITY MODELS 

A complete description of the failure models used in CARES was given by Nemeth et al. (1989) 
and (1990). These models which describe the statistical nature of fracture in brittle materials are based 
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„„ the weakest-link theory and the twtvparameter Weiboll distributor (We, bull 1939a). Concepts 
developed by Batdorf (Batdorf and Crose, 1974) provide the physical foundation for use of various failure 
criteria by incorporating linear elastic fracture mechanics into weakest-link theory to account for the 
effects of multiaxial stress states. 

The WELFEP study examines the principle of independent action ^PfA) (Barnett e t a l .1967; and 
Freudenthal, 1968), and the Weibull normal stress averaging method (Weibull, 1939b). For the P1A 
model, the probability of failure for volume flaws, Pfy, is expressed as 


Pf=l-e*p 


[-—/ 

(o7* + o^*+o^ T ) dV 


°ov 



( 1 ) 


and for surface flaws, P fg 


P*. = lexp 




( 2 ) 


where m is the Weibull modulus and a Q is the scale parameter, whose units are stress x length or 
<, a nd stress x length 2 /™ for * og . The subscripts v and s denote volume and surface flaw analysis, 
respectively. The principal stresses are represented by e v o 2 and <r 3 , V is the volume and A is e 
surface area. 

In the normal stress averaging method (shear insensitive), the probability of failure for volume flaw 
analysis is given by 


and for surface flaws, 


Pft = l -eX P **] 


(3) 


(4) 


where k is the normalized poly axial Weibull crack density coefficient and o n is the average normal 

stress. The normalized polyaxial Weibull crack density coefficient is obtained by applying the normal 
stress averaging method to a uniaxial stress state, and by equating that result with the uniaxial two- 
parameter Weibull equation. For volume flaws, k, is given by 


and for surface flaws 
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The normal stress is averaged on the surface of a unit sphere for volume flaws. Analogously, the line 
integration is performed over the contour of a unit radius circle for surface flaws. The symmetry of the 
stress reduces the integration to an eighth of the sphere and a quarter of the circle (Nemeth et al., 1990). 

The shear-sensitive models available in CARES are based on concepts developed by Batdorf and 
Crose (1974). They require a user-selected flaw geometry and a mixed-mode fracture criterion. Flaws 
intersecting the surface can be modeled as Griffith cracks, Griffith notches, or semicircular cracks. Imper- 
fections embedded in the volume can be described as Griffith flaws, Griffith cracks, or penny-shaped 
cracks. The total strain energy release rate theory is used as a mixed-mode fracture criterion for coplanar 
crack extension (Batdorf and Heinisch, 1978). Out-of-plane crack extension criteria are approximated by 
an equation with a semiempirical constant (Palaniswamy and Knauss, 1978; Shetty, 1987) that models 
various crack configurations and fracture theories such as the maximum tangential stress (Erdogan and 
Sih, 1963), the maximum strain energy density criterion (Sih, 1974), and the maximum strain energy 
release rate criterion (Hellen and Blackburn, 1975). 


PROGRAM FEATURES 

Based on the previous version of CARES (Nemeth et al., 1990), which was coupled to MSC/ 
NASTRAN and ANSYS stress analyses, a next generation algorithm has been developed that allows 
access to additional finite element programs through the use of a formatted neutral data base. In the 
following section, the new features of CARES will be highlighted in comparison to the previous version. 
The original algorithm has been reorganized into the following separate modules: 

1. Finite element data interface 

2. Reliability evaluation algorithm. 

The interface program interprets the data necessary for subsequent reliability analysis, and translates this 
information into a formatted (ASCII) neutral data base. This feature enables easy data transfer between 
different computer systems. The reliability evaluation algorithm performs the computations necessary in 
determining structural reliability. 

Furthermore, a subelement technique has been implemented. The reliability analysis is now per- 
formed at each Gauss integration point instead of using element averaging techniques. In the context of 
finite element analysis, stresses are determined at the Gaussian integration points where the local stiffness 
matrix has been evaluated. The subelement technique implies that each Gaussian integration point 
corresponds to a subelement. The subelement volume is defined as the contribution of the integration 
point to the element volume in the course of the numerical integration procedure. 

The location of the Gaussian integration point in the natural space of the finite element, as well as 
corresponding weight functions have to be considered when the subelement volume is calculated. For this 
reason, the number of subelements in each element depends on the integration order chosen, and as a con- 
sequence, the element type. The increased number of points for stress evaluation accounts for the varia- 
tion m stress over the element. Thus, considerable improvements in the accuracy of the reliability 
analysis have been realized (the reader is referred to the numerical examples that follow). 
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Due to this subelement concept a new global structure of the algorithm ,has evolved. The modular 
structure of the program readily allows future changes and additions without undertaking large program- 
ming efforts This version of CARES also allows variations in temperature throughout the structure. 
CARES 3.0 is programmed in standard FORTRAN??. This allows easy implementation on numerous 
computer systems. The introduction of one large blank common array (instead of several fixed-dimension 
arrays) in each routine allows the user to easily adapt the program to the problem size as well « to the 
available computational resources. The adaption of the size of the blank common array can readily be 
adjusted by resetting one parameter at the beginning of the program (MTOT). In contrast to pre 
versions of CARES the code can now also be implemented on PC systems. The reliability evaluation 
algorithm requires two input files, a control file and a neutral data base. The control file defines the 
number of different materials, the temperature dependent Weibull parameters and ^gs for prmting and 
postprocessing (e.g., writing the element risk of rupture intensities to a PATRAN data file). These 
parameters are stored in the blank common array. 


The finite element results and data required for the reliability analysis are supplied by the neu ra 
data base, including the stress tensors and temperatures at the integration points. The structure of e 
neutral data base is optimized with respect to memory. The finite element data is arranged within e 
neutral data base using the following hierarchy, element groups, elements and subelements. The element 
group data contain information regarding the number of elements within the group. In addition to the 
element types available in the previous version of CARES (three-dimensional brick elements for volume 
flaw analysis, membrane elements for surface flaw analysis), two-dimensional finite elements are also 
implemented. By assuming plane stress conditions or exploiting symmetry of the structure, the size o 
the finite element model can be reduced using a two-dimensional mesh. Thus, the number of degrees o 
freedom could be substantially reduced. 


The two-dimensional modeling options available (in plane loading versus bending) depend on the 
method used to determine the volume. In the case of MSC/NASTRAN, the volume is calculated by 
multiplying the area of a subelement (see NEUTRAL DATA BASE CONCEPT) by the thickness of the 
element. Therefore, bending effects are not taken into account because the stress is assumed to be 
constant through the thickness. 


Information pertaining to the elements include the number of subelements and the material identi- 
fication number. Finally, the subelement group data contains information regarding subvolume, stress 
tensor (i.e., the stress tensor at the integration point), and subelement temperature. The stresses and 
temperatures are assumed to be constant over the subelement considered. 

Depending on the size of the finite element mesh (number of element groups, elements and sub- 
elements in the finite element model) pointers are calculated for the blank common array, and the 
corresponding size limits are checked. As the stresses are transformed to principal stresses immediately 
after reading the stress tensor at the subelement level, the required memory is optimized. Only the 
principal stresses in each subelement, the subelement volume (area in the case of shell elements), the 
subelement temperature, the material identification number for the considered element and the element 
thickness (only for shell subelements) have to be stored. 

To take into account variations in the temperature field, the Weibull parameters are calculated by 
linear interpolation between the temperature dependent parameters entered in the control file. In 
CARES2, the temperature dependent material parameters were calculated by Lagrangian polynomials, 
but this algorithm did not work properly for special data configurations. 
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After performing the reliability analysis, the element risk of rupture dat« (a local measure of relia- 
bility) can be written into a PATRAN element file. This procedure allows visualization of critical regions 
in the structure (only available, if the finite element mesh has been generated via PATRAN preprocess- 
ing). For a subelement over which the state of stress and temperature are assumed to be constant, the 
reliability is expressed as 


Ri = exp [-<Pi V 4 ] 


( 7 ) 


where Rj is the reliability, <p 4 is the risk of rupture intensity and V. is the volume of the i’th sub- 
element. The reliability of the component is determined by the product of all individual subelement 
reliabilities or 


k 


r = n *i = e *p 

i«l 


v “ 


i*l 


( 8 ) 


The risk of rupture intensity for an element, (p^, is a weighted average of the intensities of the sub- 
elements or 


D 



( 9 ) 


where n is the number of subelements within the element. For shell elements (surface flaws) the volume 
is replaced by an area term. The risk of rupture intensity, 9^, is only used for PATRAN 
postprocessing. 


NEUTRAL DATA BASE CONCEPT 

The concept of a formatted neutral data base allows for the interfacing of CARES 3.0 to several finite 
element packages. The results from the finite element analysis (available in the form of standard output 
files and/or plot data files) have to be interpreted by an interface program. Due to the different formats 
used for output of the results in each finite element program this interpreter program has to be adapted 
to the finite element software used. The interface programs for NASTRAN and ABAQUS have already 
been implemented, interfaces for MARC and ANSYS are being prepared. 

The finite element results required for the reliability estimation are assembled into the formatted 
neutral data base. A further advantage of this approach is the ease of transfer of this formatted neutral 
database (ASCII file) to different computer systems. Interpreter control parameters, such as the number 
of element groups to be considered in the reliability analysis, flags for accounting for temperature field 
gradients, etc., are usually input via an interpreter control file. 

The introduction of comment lines (defined by “COM” as the first three characters of an input line) 
allows for a more friendly substructuring of the neutral data base. Based on the global hierarchy outlined 
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above, the neutral file data is arranged into records with standard FORTRAN formats: the first record 
(format A80) specifies the title for the particular reliability evaluation. The next entry (format 215) con- 
tains the number of element groups NUMEL (both volume as well as surface (shell) elements) to be con- 
sidered in the reliability analysis and the flag NUMELB indicating combined volume as well as surface 
flaw reliability analysis for shell elements. 

Depending on the parameter NUMEL, a loop over all element groups is started with the next state- 
ment (format 415): the identification number of the element group considered (INUEG), the element type 
(IEGTYP) corresponding to PATRAN element type codes, the number of elements within this element 
group (IGREL) and a shell element type indicator (ISHELL) have to be input. If shell elements are 
used in the analysis, the indicator ISHELL has to be set to 1. Information pertaining to the element 
level is read in a loop over the number of elements (IGREL) within the element group considered (format 

415, 3F 15.0). 

The element level data contains the identification number of the element considered (IELNUM), the 
number of subelements within this element (ISUBEL), the material identification number for volume 
flaw analysis (MATINP), the material identification number for surface flaw analysis (MAT2), the ele- 
ment volume (ELVOL), the element thickness (ELTHIC) only in case of shell elements, and the aver- 
aged temperature of this element (ELTEMP). The identification number of the element should be the 
same as the corresponding element number in the finite element mesh, if further postprocessing of the 
reliability results are planned. The number of subelements within this element depends on the Gaussian 
integration order chosen in the finite element analysis and on the element type. If stress data is only 
available at the center of the element (standard case in NASTRAN analyses prepared for previous 
CARES versions), one subelement is specified (ISUBEL = 1). The material identification numbers must 
correspond to the material numbers used in the finite element analysis as well as the material numbers on 
the control file (containing the Weibull parameters). Element volume and element temperature are 
necessary only for element group information summaries in the printout file. For the reliability evalu- 
ation the subelement data are directly used. 

A third loop is started for reading the subelement level data. In case of volume type elements 
(ISHELL = 0 on Element Group Entry) the identification number of the subelement (ISUNUM), the 
subelement volume (SUBVOL), the subelement temperature (SUBTEM) and the full stress tensor at the 
integration point (defined by o , <r xy , o yi and in the local coordinate system) have to be 

specified (format I5,5F15.0,/,35X,3F15.0). 

In the case of shell type elements (ISHELL = 1 on Element Group Entry) plane stress conditions are 
assumed. For that reason only o and a xy are necessary for the definition of the subelement 

stress tensor. Instead of the subelement volume, the subelement area has to be specified (15, 5F 15.0): 
ISUNUM, SUBAREA, SUBTEM, a ^ and o^. 

Typically, subelement volumes are not included with standard finite element output. Thus, the vol- 
ume of each subelement (corresponding to a Gauss integration point) is calculated in the interpreter pro- 
gram using the shape functions inherent to the element type. In the usual context of finite element 
methods, the volume of a three-dimensional element, (i.e., brick, wedge, pyramid or tetrahedron) is calcu- 
lated after transformation into the natural coordinate space (Bathe, 1982), 

1 1 1 

V= f f f det J(r^) drdsdt ( 10 ) 

-l -l -l 
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where J is the Jacobian operator and r,s,t are the natural coordinates. By applying a Gaussian inte- 
gration scheme, the integral can be expressed as 


v= E £ £ ** J (w*) w, Wj w k (ii) 

i-1 j-1 k-1 

where n is the integration order, rj,Sj,t k are the coordinates of the Gaussian integration points in 
natural space and the W’s are the weight functions. The area of an element (i.e., quadrilateral or tri- 
angle) is determined in a similar manner 

A= £ £ ** w , w, ( 12 ) 

i-l j-1 


The volume of an axisymmetric element can also be determined similarly 

v= R £ £ ** J (vo w, w. 

i-i j-i 

where R is the radial coordinate of the integration point in the global coordinate system. 


(13) 


NUMERICAL EXAMPLES 

To demonstrate the improvements in the algorithm, example problems are selected from a round- 
robin conducted by the WELFEP group. The first test case is a beam in pure bending. Comparisons 
with the analytical solution are made to verify convergence and the mathematical correctness of the 
codes. The second structure is a notched beam in four-point bending. The induced concentration will 
test the algorithms used to identify element stresses. The final test case is a tube under ideal torsion. 
Comparisons from this problem verify the manner in which the reliability for multiaxial stress states is 
determined. For the purpose of this study, two failure models will be examined: the Principle of Inde- 
pendent Action (PIA) and the Normal Stress Averaging method (NSA). Both are based on WeibulPs 
weakest link theory. 

To demonstrate the improvements achieved in this version of CARES 3.0, the results of the reliability 
analysis are compared to the reliability data generated by the previous version of the program (CARES2: 
subelements are obtained by dividing the element into equal cells, the stresses are interpolated from the 
nodes.). MSC/NASTRAN finite element analysis is used for each of the test cases. CARES 3.0 is used 
to determine the reliability for volume and surface flaw analysis. All volumes are meshed with three- 
dimensional elements (CHEXA) and shell elements (CQUAD8) are used for modeling the surfaces (mid- 
side nodes are present). Material properties are arbitrarily chosen (the same for each test case: Young’s 
modulus is 300 MPa and Poisson’s ratio is 0.2). 

A description of the beam analyzed in test case number 1 is given in Fig. 1. Bending of the beam is 
induced by prescribed displacements that are linearly varied at x = L. Due to symmetry considerations 
only one quarter of the beam has to be modeled (as indicated in Fig. 1). In order to study the influence 
of finite element discretization, linear finite element analyses are performed with a coarse and a fine mesh. 
The coarse mesh consists of 10 elements along the length direction (x), four in the height direction (y), 
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and one element in the thickness direction (z). The fine mesh is similar with 20, 8 and one element in the 
length, height, and thickness directions, respectively. The results of the CARES 3.0 reliability analyses 
as well as the analytical solution are listed in Table I. 

As expected, the solutions for the refined mesh are much closer to the analytical value than the 
results for the coarse mesh. Comparing the results from CARES 3.0 to the results for the previous ver- 
sion CARES2 shows significant improvements in accuracy. Even in the case of the coarse mesh, the 
increased number of subelements enables the model to approach the exact solution. The results for the 
fine mesh nearly coincide with the analytical solution. By modeling the component with a coarse mesh, 
numerical efforts in finite element analysis can be reduced for preliminary design. The final design is 
usually checked with a fine mesh. 

In Example 2, a notched beam in four-point bending is investigated. Due to the notch, a stress con- 
centration is induced which yields a multidimensional stress field. Symmetry considerations again allow 
us to reduce the finite element model to one quarter of the original beam (as shown in Fig. 2). In the 
vicinity of the notch the finite element model is refined to account for the strongly varying stress gradi- 
ent. The reliability results are listed in Table I. An analytical solution for the reliability is not available. 
Comparing the results from CARES 3.0 with the results from CARES2 shows that the new version calcu- 
lates more conservative results for volume flaws, and lower values for probability of failure for surface 
flaws. 

In Example 3, a tube under pure torsion is analyzed. Figure 3 shows the geometry of the tube. Tor- 
sional stresses are induced by prescribed displacements at the end of the tube (z = L). Three elements 
are used in the radial direction, 24 around the circumference, and only one element in the length direc- 
tion. The results of the reliability analysis are shown in Table I, including the analytical solution of this 
problem for the PIA model. Again, CARES 3.0 delivers more conservative results than CARES2. The 
probability of failure data coincides with the analytical solution for the PIA model. To achieve the same 
degree of accuracy with CARES2, a refined mesh would be necessary. Thus, CARES 3.0 allows con- 
siderable reduction in the size of finite element discretizations. 


SUMMARY 

Based on the previous version, CARES2, an improved analytical tool has been developed for the reli- 
ability evaluation of monolithic ceramic components. By dividing the original program into two separate 
modules, a greater flexibility is achieved leading to easy integration with various finite element programs. 
Additionally, the neutral data base concept allows easy data transfer between different computer systems. 
The modular structure also allows for the incorporation of more sophisticated failure models. CARES 3.0 
will serve as the platform for future time dependent analysis. 

With the introduction of a new subelement technique, the accuracy of the reliability analyses has 
improved. In order to achieve the same degree of accuracy, CARES 3.0 allows structural components to 
be modeled with coarser finite element meshes than what was required by the previous version of the pro- 
gram. Numerical examples demonstrate this feature. 
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APPENDIX A 


PREPARATION OF A MSC/NASTRAN PROBLEM FOR RELIABILITY ANALYSIS 

WITH CARES 3.0 

The following describes the necessary information to be included in the MSC/NASTRAN input for 
further analysis by CARES 3.0. Figure Al shows the operational flow for the reliability evaluation of 
ceramic components. This analysis includes the CARES integrated design program, MSC/NASTRAN 
static finite element analysis and PATRAN PLUS pre- and post-processing. NASCARES interprets 
MSC/NASTRAN output and assembles it into the neutral file. 

The MSC/NASTRAN data file contains three parts: the EXECUTIVE CONTROL, CASE CON- 
TROL and BULK DATA. The following example of the EXECUTIVE CONTROL is typical for a 
problem involving static analysis; 

NASTRAN DAYLIMIT = -1 

ID CERAMIC, FRACTURE 

SOL 61 

APP DISP 

TIME 30 

CEND 

It is assumed that the user is familiar with MSC/NASTRAN, so each statement will not be explained 
herein. For problems that are modeled with cyclic symmetry, Solution Sequence 61 is replaced by Rigid 
Format 47. 

Previously, Solution Sequence 61 (SUPERELEMENT STATICS or in version 66, SOL 101) and 
Rigid Format 47 (CYCLIC STATICS) were used for the analysis. These solution methods were 
employed because element volumes and areas can be obtained through a PARAM card located in the 
BULK DATA. Since the element (and subelement) volumes/areas are now calculated by NASCARES, 
there is no longer a limit on solution sequences selection. 

For the CASE CONTROL DECK, a typical problem may contain, 

TITLE = TITLE 

SUBTITLE = SUBTITLE 

SEALL = ALL 

SPC=50 

MPC=51 

LOAD=20 

STRESS(PRINT)=ALL 

ECHO=SORT 


BEGIN BULK 

The following three items are important to note. First, multiple load subcases are not defined. CARES 
3.0 does not handle stress output from multiple subcases. Therefore, separate NASTRAN executions are 
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required for each loading condition. Secondly, element stresses are routed to the printout file when 
STRESS(PRINT)=ALL; is specified. All element stresses must be placed in the print file. Lastly, 
ECHO=SORT; will send the sorted BULK DATA to the print file. NASCARES reads the necessary 
BULK DATA, such as material ID, element connectivities, and nodal temperatures, from the print file. 
The BULK DATA entries which are processed by NASCARES are listed in table Al. 

The PARAM,EST,1 card will output element volumes and areas to the NASTRAN printout file. 

This card only functions for some solution sequences or rigid formats. The TEMPD card is required if 
some grid points do not have an explicit temperature assignment. All nodes must have an associated 
temperature. These are specified with TEMP, TEMP*, or TEMPD card. Since this version of CARES 
does not examine cards for their subcase identity, caution must be exercised not to allow multiple temper- 
ature cases. Nodal temperatures are used by CARES to obtain the average element or subelement tem- 
peratures, which are not routinely available in MSC/NASTRAN. The CARES program utilizes results 
from isoparametric two- and three-dimensional finite elements. A list of the available elements is given in 
table A2. For typical designs, the three-dimensional HEXA and PENTA elements should be used to 
model the structure, including thin cross sections such as blades. 

After execution of the MSC/NASTRAN problem, the analyst should have a print file containing the 
sorted BULK DATA, element stresses and general problem information. At this point, all the informa- 
tion required from MSC/NASTRAN to determine the component reliability is present. 
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APPENDIX B 


NEUTRAL FILE SETUP 


PURE BENDING BY PRESCRIBED DISPLACEMENTS 
COM 

COM EXAMPLE NEUTRAL FILE 
COM 


2 0 

COM FORMAT (215) 

COM NUMEL = 2 : NUMBER OF ELEMENT GROUPS 

COM NUMELB = 0 : FLAG FOR SHELL ELEMENT GROUPS WITH VOLUME PROCESSING 

COM 0 = NO VOLUME PROCESSING OF SHELL ELEMENT GROUPS 

COM 1 = VOLUME PROCESSING OF SHELL ELEMENT GROUPS 

COM 


COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 

COM 


COM 

COM 

COM 

COM 

COM 

COM 

COM 


1 8 40 0 

FORMAT (415) 

INUEG = 1 : GROUP NUMBER 

IEGTYP = 8 : ELEMENT TYPE — EQUIVALENT TO PATRAN’S SHAPE PARAMETER 

3 = TRIANGLE 

4 = QUADRILATERAL 

5 = TETRAHEDRON 

6 = PYRAMID 

7 = PENTAHEDRON (WEDGE) 

8 = HEXAHEDRON (BRICK) 

IGREL = 40 : NUMBER OF ELEMENTS IN THIS GROUP 
ISHELL = 0 : SHELL ELEMENT FLAG 

0 = VOLUME OR AXISYMMETRIC ELEMENT GROUP 

1 = SHELL ELEMENT GROUP 


27 300 


FORMAT 
IELNUM = 
ISUBEL = 
MATINP = 
MAT2 = 
ELVOL = 
ELTEMP = 
ELTHIC = 


(4I5,3E15.6) 


0 . 500000E+00 0.700000E+02 O.OOOOOOE+OO 


1 
27 
300 
0 

.50000 

70.0 

0.0 


ELEMENT NUMBER 

NUMBER OF SUBELEMENTS 

MATERIAL NUMBER FOR VOLUME FLAW ANALYSIS 

MATERIAL NUMBER FOR SURFACE FLAW ANALYSIS 

ELEMENT VOLUME 

ELEMENT TEMPERATURE 

ELEMENT THICKNESS (SET TO ZERO FOR VOLUME ELEMENTS) 


0. 107167E-01 0 . 700000E+02 -0. 166906E+03 

0.201770E-06 

FORMAT (15 , 5E15 . 6 , / , 35X , 3E15 . 6) 


ISUNUM = 1 
SUBVOL = .01072 
SUBTEM = 70.0 
SSTR(L) ,L=1,6 


SUBELEMENT NUMBER 
SUBELEMENT VOLUME 
SUBELEMENT TEMPERATURE 
STRESSES a 


xy 


0 . 144637E-05 -0 . 1 67879E-06 

0 . 137424E-05 -0 . 204341E-06 




2 0. 171468E-01 

3 0. 107167E-01 

4 0.171468E-01 

5 0.274348E-01 


0.700000E+02 
0.700000E+02 
0 . 700000E+02 
0.700000E+02 


-0. 166906E+03 
0.719933E-06 
-0. 166906E+03 
0.123810E-05 
-0 . 225001E+03 
-0. 113103E-07 
-0 . 225001E+03 


-0.128449E-05 
-0. 108985E-05 
-0.401535E-05 
-0.355393E-05 
0 . 500609E-06 
0.451752E-06 
-0.753135E-06 


-0.732417E-07 
-0.536754E-06 
0.213957E-07 
-0.869168E-06 
-0. 196162E-06 
-0. 104590E-06 
-0. 166012E-06 


12 



0 . 179052E-07 -0 . 2936840-06 -0 . 277750E-06 


25 

0.107167E-01 

0 . 700000E+02 

26 

0.171468E-01 

0 . 700000E+02 

27 

0.107167E-01 

0.700000E+02 


0 . 1 69347E+02 0 . 780221E-01 

-0 . 536298E+00 -0 . 645437E+00 
0 . 173446E+02 0 . 174948E+01 

0 . 330185E+00 0 . 908340E+00 

0 . 177545E+02 0 . 342093E+01 

0. 119667E+01 0.246212E+01 


0.234190E-01 
0.990120E-01 
0 . 429904E+00 
0 . 978093E-01 
0 . 836389E+00 
0 . 966066E-01 


4 64 1 

FORMAT (415) 


3TYP = 4 : ELEMENT TYPE — EQUIVALENT TO PATRAN’S SHAPE PARAMETER 

JEL = 64 : NUMBER OF ELEMENTS IN THIS GROUP 
JELL = 1 : SHELL ELEMENT FLAG 

0 = VOLUME OR AXISYMMETRIC ELEMENT GROUP 

1 = SHELL ELEMENT GROUP 

9 0 300 0 . 500000E+00 0.700000E+02 0.100000E-05 

FORMAT (4I5,3E15.6) 

LNUM = 2 : ELEMENT NUMBER 

JBEL = 9 : NUMBER OF SUBELEMENTS 

riNP = 0 : MATERIAL NUMBER FOR VOLUME FLAW ANALYSIS 

12 = 300 : MATERIAL NUMBER FOR SURFACE FLAW ANALYSIS 

VOL = .50000 : ELEMENT AREA 

TEMP = 70.0 : ELEMENT TEMPERATURE 

THIC = l.OE-6 : ELEMENT THICKNESS 

0.385802E-01 0.700000E+02 -0 . 300001E+03 -0.125358E-06 0.202518E-06 

FORMAT (I5,5E15.6) 

UNUM = 1 : SUBELEMENT NUMBER 

BVOL = .03858 : SUBELEMENT AREA 

BTEM = 70.0 : SUBELEMENT TEMPERATURE 

TRfLI .L=l .3 : STRESSES o r 


2 0.617284E-01 

3 0.385803E-01 

4 0.617284E-01 

5 0.987654E-01 

6 0.617284E-01 

7 0.385803E-01 

8 0.617284E-01 

9 0.385803E-01 

3 9 0 300 

1 0.771605E-01 

2 0.123457E+00 


0 . 700000E+02 
0 . 700000E+02 
0 . 700000E+02 
0 . 700000E+02 
0 . 700000E+02 
0.700000E+02 
0.700000E+02 
0 . 700000E+02 
0. 100000E+01 
0 . 700000E+02 
0.700000E+02 


-0 . 300001E+03 
-0 . 300001E+03 
-0.300001E+03 
-0 . 300001E+03 
-0 . 300001E+03 
-0.300001E+03 
-0.300001E+03 
-0.300001E+03 
0.700000E+02 
-0.166906E+03 
-0. 166906E+03 


-0.341787E-06 
-0. 55821 6E-06 
-0 . 335760E-06 
0.868193E-07 
0 . 509398E-06 
-0.546161E-06 
0.515426E-06 
0.157701E-05 
0.100000E-05 
-0.823890E-06 
0.719767E-06 


0.765522E-06 
0 . 132853E-05 
0.251547E-06 
0.837670E-06 
0.142379E-05 
0.300575E-06 
0.909818E-06 
0. 151906E-05 

0 . 366659E-06 
-0.230519E-06 


5 0.987654E-01 

6 0.617284E-01 

7 0.385802E-01 

8 0.617284E-01 

9 0.385803E-01 


0 . 700000E+02 
0.700000E+02 
0 . 700000E+02 
0 . 700000E+02 
0.700000E+02 


-0. 183567E+02 
-0. 181104E+02 
-0.438719E+01 
-0.381854E+01 
-0.324989E+01 


-0 . 186999E+02 
-0.171232E+02 
-0.710970E+01 
-0.418869E+01 
-0.126768E+01 


0.887556E+00 
0. 176388E+01 
0.386692E-01 
0. 195399E+01 
0.386931E+01 
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TABLE 1— FAILURE PROBABILITIES FOR THE 1 


WELFEP TEST CASES 


Fracture 

criteria 

Mesh 

E I 


CARES 2 

CARES 3.0 

Analytical 

solution 

1 — Pure bending by prescribed displacement 

s 

PIA 

Volume 

Coarse 

Fine 

Eh 


0.1364 

Surface 

Coarse 

Fine 

n 

0.7677 

.7691 

0.7691 

NSA 

Volume 

Coarse 

Fine 

0.1053 

.1278 

0.1336 

.1363 

0.1364 

Surface 

Coarse 

Fine 

0.7524 

.7647 

0.7679 

.7693 

0.7691 

2 — Notched beam in four-point bend 

PIA 

Volume 

(a) 

0.1248x10“* 

0.1304xl0“ 4 


Surface 

(a) 

0.1184x10“* 

0.1043x10“* 


NSA 

Volume 

(a) 

0.1340x10“ 4 

0.1421xl0~ 4 


Surface 

M 

0.1209X10" 3 

0.1065x10“* 


3 — Pure torsion by prescribed displacements 

PIA 

Volume 

(a) 

0.1892 

0.1933 

0.1933 

Surface 

(a) 

0.4891 

0.4891 

0.4891 

NSA 

Volume 

(a) 

0.1371 

0.1402 


Surface 

(a) 

0.3757 

0.3757 



% Mesh u defined by WELFEP proposa.li. 
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TABLE Al. — KEY 
BULK DATA 


CHEXA 

C0RD2C 

C0RD2R 

C0RD2S 

CPENTA 

CQUAD4 

CQUAD8 

CTETRA 

CTRIA3 

CTRIA6 

CTRIAX6 

GRDSET 

GRID 

PSHELL 

PSOLID 

TEMP 

TEMPD 


TABLE A2. — MSC/NASTRAN SUPPORTED 
ELEMENT TYPES 


(Shell elements should be defined to have membrane proper- 
ties (i.e., no bending or transverse shear) and be very thin 
(i.e., l.E-06) so as not to contribute additional stiffness.] 



Element 

type 

Number 
of nodes 

Number of 
subelements 
(if divided) 

Solid 

HEXA 

8-19 

8 



20 

27 


PENTA 

mm 

6 



■EH 

9 


TETRA 

4-9 

1 



10 

4 

Shell 

QUAD4 

4 

1 


QUAD8 

mm 

4 




9 


TRIA3 

3 

1 


TRIA6 

3-5 

1 



6 

3 

Axisymmetric 

TRIAX6 

3-5 

1 



6 

3 
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Figure A1 - Block diagram for the analysis and fast-fracture reliability evaluation 
of ceramic components. This analysis includes the CARES integrated 
design program, MSC/NASTRAN finite element analysis and PATRAN PLUS pre- 
and post-processing. 
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Data 

L = 10.0 mm 
H = 4.0 mm 
B = 1.0 mm 
m = 15.0 

o ov = 362.4 MPa x mm 0 - 200 
o os = 362.4 MPa x mm 0 - 133 

Figure 1 .—Pure bending by prescribed displacements. 
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Data 

L, = 10.0 mm 
l_ 2 = 20.0 mm 
L 3 = 22.5 mm 
H * 4.0 mm 
B s 1.0 mm 
m = 15.0 

o ov = 1035.6 MPa x mm 0 - 200 
a os = 1035.6 MPa x mm 0 - 133 

L N = 1 .0mm 
R n = 0.5mm 

Figure 2.— Notched beam in 4-pt bend. 
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Data 

L * 12.5 mm 
Dj = 1 0.0 mm 
D 0 = 12.0 mm 
m = 15.0 

a w «1 55.3 MPa x mm 0 - 200 
Oos - 1 55.3 MPa x mm 0 - 133 

Figure 3. — Pure torsion by prescribed displacements. 
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